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We present a study of the D + ir~ , D°ty + , and D* + -k~ systems in inclusive e + e~ — > cc interactions 
in a search for new excited D meson states. We use a dataset, consisting of ~454 fb _1 , collected at 
center-of-mass energies near 10.58 GeV by the BABAR detector at the SLAC PEP-II asymmetric- 
energy collider. We observe, for the first time, candidates for the radial excitations of the D°, D*°, 
and D* + , as well as the L = 2 excited states of the D° and D + , where L is the orbital angular 
momentum of the quarks. 

PACS numbers: 14.40.Lb, 13. 25. Ft, 12.38.-t 



The spectrum of mesons consisting of a charm and an 
up or a down quark is poorly known. The spectrum of 
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quark- ant iquark systems was predicted in 1985 using a 
relativistic chromodynamic potential model [lj . The low- 
mass spectrum of the cu or cd system is comprised of the 
ground states (IS), the orbital excitations with angular 
momentum L=l,2 (1P,1D), and the first radial excita- 
tions (2S). In this paper we label the states using the no- 
tation D^j S+1 \nL), where J is the total angular momen- 
tum of the state, n is the radial quantum number, and L 
and S are the orbital angular momentum and total spin 
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of the quarks. Besides the ground states (D,D*), only 
two IP states, known as the D x (2420) and #2(2460) 0, 
are well-established experimentally since they have rela- 
tively narrow widths (~30 MeV). In contrast, the other 
two IP states, known as the .0(5(2400) and £>i(2430), are 
very broad (^300 MeV), making them difficult to detect 



To search for states not yet observed, we analyze the 
inclusive production of the Z) + 7r~, D°ir + , and D* + n~ [6\ 
final states in the reaction e + e~ — > cc — !> D^ttX, where 
X is any additional system. We use an event sample con- 
sisting of approximately 590 million e + e~ — > cc events 
(454 fb" 1 ) produced at e + e~ center-of-mass (CM) ener- 
gies near 10.58 GeV and collected with the BaBAR de- 
tector at the SLAC PEP-II asymmetric-energy collider. 
Our signal yield for the L — 1 resonances is more than ten 
times larger than the best previous study Q , resulting in 
much greater sensitivity to higher resonances. 

The BaBAR detector is described in detail in Ref. Q. 
Charged-particle momenta are measured with a 5-layer, 
double-sided silicon vertex tracker (SVT) and a 40- 
layer drift chamber (DCH) inside a 1.5-T supercon- 
ducting solenoidal magnet. A calorimeter consisting of 
6580 CsI(Tl) crystals is used to measure electromagnetic 
energy. A ring-imaging Cherenkov radiation detector 
(DIRC), aided by measurements of ionization energy loss, 
dE/dx, in the SVT and DCH, is used for particle identi- 
fication (PID) of charged hadrons. 

The Dir system is reconstructed in the D + ir~ and 
D°ir + modes, where D + ->• K~tt + tt + and D° ->• K~n + . 
A PID algorithm is applied to all tracks. Charged kaon 
identification has an average efficiency of 90% within the 
acceptance of the detector and an average pion-to-kaon 
misidentification probability of 1.5%. 

For all channels we perform a vertex fit for the D + 
and D° daughters. To improve the signal to background 
ratio for D + — > K~Tt + n + , we require that the measured 
flight distance of the D + candidate from the e + e~ in- 
teraction region be greater than 5 times its uncertainty. 
To improve the signal purity for D° — > K~ir + we re- 
quire cos 8k > — 0.9 where 9 k is the angle formed by the 
K~ in the D° candidate rest frame with respect to the 
prior direction of the D° candidate in the CM reference 
frame. The Dir candidates for both D + and D° are then 
reconstructed by performing a vertex fit with an addi- 
tional charged primary pion, which originates from the 
e + e~ interaction region. For all vertex fits we require a 
X 2 probability > 0.1%. 

In the D°ir + sample, we veto D° candidates from D* + 
or D*° decays by forming D°tt + (where the ir + is any 
additional pion in the event) and D°ir° combinations, 
and rejecting the event if the invariant-mass difference 
between this combination and the D° candidate is within 
2(7 of the nominal D*-D mass difference [2|, where a is 
the detector resolution. 

The K~Tr + n + and K~Tt + mass distributions are shown 
in Figs. Q] a) and [I] b). We fit these distributions to 
a linear background and a Gaussian signal; the signal 
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FIG. 1: (color online) Mass distribution for a) D + and b) 
D° candidates in the D + -k~ and D°-k + samples. Plots c) 
and d) correspond to the D* + ir~ sample and show the mass 
distribution for D° candidates and the Am distribution for 
D* + candidates. The vertical lines show the signal and, in a) 
and b), the side- band regions. 



widths obtained are <j D + — 6.7 MeV/c 2 and <tjjo = 7.6 
MeV/c 2 . The signal region is defined to be within ±2.5ct 
of the peak, while sideband regions are defined as the 
ranges (±5.0cr, ±7.5cr) and (±4.0cr, ±6.5cr) for the D + 
and D°, respectively. The D + signal region has purity 
N S /{N S + N B ) = 65%, where N s (N B ) is the number of 
signal (background) events, while the D° purity is 83%. 

The D* + ir~ system is reconstructed using the D a — > 
K~tt + and D° — > K^tt + tt^tt + decay modes. A D° 
candidate is accepted if its invariant mass is within 30 
MeV/c 2 of the mean value. A D* + candidate is re- 
constructed by requiring an additional slow pion (7r+) 
originating from the e + e _ interaction region. We se- 
lect a D* + candidate if the mass difference Am = 
m{K~ tt + (n + tt~ )nf ) — m(K~ir + (n + ir~)) is within 2.0 
MeV/c 2 of the mean value. The D° candidate invariant 
mass distribution and the Am distribution are shown in 
Figs. CD c) and [T]d). The D*+ signal purity is 89%. Fi- 
nally, we reconstruct a D* +/ n~ candidate by combining 
a D* + candidate with an additional charged track iden- 
tified as a 7r~ and applying a vertex fit. 

Background from e + e~ — > BB events, and much of 
the combinatorial background, are removed by requiring 
the CM momentum of the D^tt system to be greater 
than 3.0 GeV/c. In addition, we remove fake primary 
pion candidates originating mainly from the opposite side 
of the event by requiring cos 9^ > —0.8. The angle 0„ 
is defined in the D^tt rest frame as the angle between 
the primary pion direction and the prior direction of the 
DWtt system in the CM frame. 

To extract the resonance parameters we define the vari- 
ables M(D + ir~) = m(K^Tr + TT + TT^) - m(if~7r + 7r + ) + 
m D + and M(L> 7r+) = m(K ~7r+7r+) - m(K~n + ) +m D a, 
where m D + and m D o are the values of the D + and D° 
mass @- The use of the mass difference improves the 
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resolution on the reconstructed mass to about 3 MeV/c 2 . 
We remove the contribution due to fake D + and D° can- 
didates by subtracting the M(Dn) distributions obtained 
by selecting events in the D + or D° candidate mass side- 
bands. 

The D + tt~ and D°ir + mass spectra are presented in 
Fig. [5] and show similar features. 

• Prominent peaks for D^(2460) and £^(2460)+. 

• The D + tt~ mass spectrum shows a peaking back- 
ground (feeddown) at about 2.3 GeV/c 2 due to de- 
cays from the £>i(2420)° and £^(2460)° to D* + Tr~. 
The D* + in these events decays to D + ir° and the 7r° 
is missing in the reconstruction. The missing ir° has 
very low momentum because the D* + decay is very 
close to threshold. Therefore, these decays have a 
mass resolution of only 5.8 MeV/c 2 and a bias of 
-143.2 MeV/c 2 . Similarly, D Q n + shows peaking 
backgrounds due to the decays of the Z?i(2420) + 



xiooo 



and DjSi(2460) H 



decays to 



D U 7T 



• Both D + tt and D°ir + mass distributions show 
new structures around 2.6 and 2.75 GeV/c 2 . We 
call these enhancements D*(2600) and £>*(2760). 

We have compared these mass spectra with those ob- 
tained from generic e + e~ —> cc Monte Carlo (MC) 
events. These events were generated using JETSET Q 
with all the known particle resonances incorporated. The 
events are then reconstructed using a detailed GEANT4 
fioj detector simulation and the event selection proce- 
dure used for the data. In addition, we study Dtt mass 
spectra from the D + and D° candidate mass sidebands, 
as well as mass spectra for wrong-sign D + ir + and D°ir~ 
samples. We find no backgrounds or reflections that can 
cause the structures at 2.6 and 2.76 GeV/c 2 . In the study 
of the D°ir + final state we find a peaking background 
due to events where the D° candidate is not a true D°, 
but the K~ candidate and the primary 7r + candidate are 
from a true D° — > K~ir + decay. These combinations 
produce enhancements in M(D°ir + ) both in the D° can- 
didate mass signal region and sidebands. However, this 
background is linear as a function of the D° candidate 
mass, is removed by the sideband subtraction. 

The smooth background is modeled using the function: 



B{x) = P{x) x 



O do+d 1 x+d 2 x 2 



for X < Xq, 
for X > Xq, 



(1) 



where P(x) = -^.\J\x 2 — (jno + m^) 2 }^ 2 — (mjj — m^) 2 } 
is a two-body phase-space factor and x = M(Dir). Only 
four parameters are free in the piece-wise exponential: 
Ci , C2 , d,2 , and xq . The parameters do and di are fixed by 
requiring that B(x) be continuous and differentiable at 
the transition point xq. We account for the feeddown of 
peaking backgrounds by convolving Breit-Wigner (BW) 
functions [11] with a function describing the resolution 
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FIG. 2: (color online) Mass distribution for D + ir~ (top) and 
D°ir + (bottom) candidates. Points correspond to data, with 
the total fit overlaid as a solid curve. Fhe dotted curves are 
the signal components. The lower solid curves correspond 
to the smooth combinatoric background and to the peaking 
backgrounds at 2.3 GeV/c 2 . The inset plots show the distri- 
butions after subtraction of the combinatoric background. 



and bias obtained from the simulation of these decays. 
The mass and width of the Z)i(2420) feeddown are fixed 



to the values obtained in the D 



analysis described 



below, while the parameters of the D^IAQQ) feeddown 
are fixed to those of the true Z?2(2460) in the same 
M(Dtt) distribution. 

The DJ(2460) is modeled using a relativistic BW func- 
tion with the appropriate Blatt-Weisskopf centrifugal 
barrier factor 0. The I?* (2600) and L>*(2760) are mod- 
eled with relativistic BW functions |2;J . Finally, although 
not visible in the M{D + n~) mass distribution, we in- 
clude a BW function to account for the known resonance 
Dq(2400), which is expected to decay to this final state. 
The x 2 P er number of degrees of freedom (NDF) of the fit 
decreases from 596/245 to 281/242 when this resonance 
is included. This resonance is very broad and is present 
together with the feeddown and Z?2(2460)°; therefore we 
restrict its mass and width parameters to be within 2a 
of the known values [5| . The shapes of the signal compo- 
nents are corrected for a small variation of the efficiency 
as a function of M(Dtt) and are multiplied by the two- 
body phase-space factor. They are also corrected for the 
mass resolution by convolving them with the resolution 
function determined from MC simulation of signal de- 
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FIG. 3: (color online) Mass distributions for D* + tt~ candi- 
dates. Top: candidates with |cos#h| > 0.75. Middle: can- 
didates with |cos#£r| < 0.5. Bottom: all candidates. Points 
correspond to data, with the total fit overlaid as a solid curve. 
The lower solid curve is the combinatoric background, and 
the dotted curves are the signal components. The inset plots 
show the distributions after subtraction of the combinatoric 
background. 



cays. The fit to the M(D + tt~) distribution (Fit A) is 
shown in Fig. [5] (top). The results of this fit, as well as 
fits to the other final states described below, are shown 
in Table |TJ In this table, the significance for each new 
signal is defined as the signal yield divided by its total 
uncertainty. 

The fit to the D°ir + mass spectrum is similar to that 
described for the D + n~ system. Because the feeddown 
is larger and the statistical precision of the resonances is 
not as good as for D + ir~, we fix the width parameters 
of all resonances to the values determined from D + tt~ 



assuming isospin symmetry. The fit to the M(D°tt + ) 
mass distribution (Fit B) is shown in Fig. [2] (bottom); 
this fit has x 2 /NDF of 278/224. We find consistent mass 
values for both D*(2600) and D*(2760) in the fits of the 
D + ir~ and D°ir + mass distributions. 

We now search for these new states in the D* + tt~ 
decay mode. We define the variable M(D* + %~) = 
m(K~ir + (ir + Tr~)irfTr~) — m{K ~ ir + (tt + tt~ )irf ) + m D ,+ 
where m D *+ is the value of the D* + mass [4]. The D* + ir~ 
mass distribution is shown in Fig. [3] and shows the fol- 
lowing features. 

• Prominent L»i(2420)° and D|(2460)° peaks. 

• Two additional enhancements at ^2.60 GeV/c 2 
and ^2.75 GeV/c 2 , which we initially denote as 
D*(2600)° and £>(2750)°. 

Studies of the generic MC simulation as well as studies of 
the D* + sidebands and the wrong-sign sample (D* + ir + ) 
show no peaking backgrounds in this mass spectrum. 

We fit M(D* + ir~) by parametrizing the background 
with the function in Eq. Q. The £> 1 (2420)° and 
£>2(2460)° resonances are modeled using relativistic BW 
functions with appropriate Blatt-Weisskopf form factors. 
The D*(2600)° and D(2750)° are modeled with relativis- 
tic BW functions. The broad resonance £^(2430)° is 
known to decay to this final state, however, this fit is in- 
sensitive to it due to its large width (~380 MeV) 0] and 
because the background parameters are free. 

Due to the vector nature of the D* + , the D* + n~ fi- 
nal state contains additional information about the spin- 
parity ( J p ) quantum numbers of the resonances. In the 
rest frame of the D* + , we define the helicity angle 9jj as 
the angle between the primary pion 7r~ and the slow pion 
7r + from the D* + decay. The distributions in cos Oh for 
the predicted resonances, assuming parity conservation, 
are given in Table HT1 Initially, we have attempted to fit 
the M(D* + tt~) distribution incorporating only two new 
signals at ~2.6 GeV/c 2 and at ~2.75 GeV/c 2 . However, 
when we extract the yields as a function of cos Oh we find 
that the mean value of the peak at ~2.6 GeV/c 2 increases 
by ^70 MeV/c 2 between cos Oh = — 1 and cos Oh = 0, 
and decreases again as cos Oh — > +1. This behaviour 
suggests two resonances with different helicity-angle dis- 
tributions are present in this mass region. To proceed we 
incorporate a new component, which we call _D(2550)°, 
into our model at ^2.55 GeV/c 2 . We extract the param- 
eters of this component by requiring |cos6>#| > 0.75 in 
order to suppress the other resonances. In this fit (Fit 
C), shown in Fig. [3] (top), we fix the parameters of the 
Z^(2460)° and L»*(2600)° to those measured in D+ij-. 
We obtain a v 2 / NDF of 214/205 for this fit. This fit also 
determines the parameters of the Di(2420)°. We then 
perform a complementary fit (Fit D), shown in Fig. 
(middle), in which we require | cos 0h\ < 0.5 to discrimi- 
nate in favor of the L>*(2600)°. We obtain a x 2 /NDF of 
210/209 for this fit. To determine the final parameters of 
the D(2750)° signal we fit the total D* + n~ sample while 
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TABLE I: Summary of the results. The first error is statistical and the second is systematic; "fixed" indicates the parameters 
were fixed to the values from Fit A or C. The significance is defined as the yield divided by its total error. 



Resonance 


Channel(Fit) 


Efficiency (%) 


Yield (xlO 3 ) 


Mass (MeV/c 2 ) 


Width (MeV) 


Significance 


Di(2420) ,J 


D* + 1T 


WJ 




102.8±1.3±2.3 


2420.1±0.1±0.8 


31.4±0.5±1.3 






D* + n~ 




f .09 ± 0.03 


214.6±1.2±6.4 


2420.1 (fixed) 


31.4(fixed) 




D5(2460)° 


D + TT~ 


(A) 


1.29 ±0.03 


242.8±1.8±3.4 


2462.2±0.1±0.8 


50.5±0.6±0.7 






D* + 7V' 


~ (E) 


1.12 ±0.04 


136±2±13 


2462.2(fixed) 


50.5(fixed) 




D(2550) u 


D* + tt' 


" (C) 




34.3±6.7±9.2 


2539.4±4.5±6.8 


130±12±13 


3.0a 




D* + k 


~ (E) 


1.14 ±0.04 


98.4±8.2±38 


2539.4(fixed) 


130(fixed) 




D*(2600)° 


D+-R- 


(A) 


1.35 ±0.05 


26.0±1.4± 6.6 


2608.7±2.4±2.5 


93±6±13 


3.9a 




D*+tv- 


" (D) 




50.2±3.0±6.7 


2608.7(fixed) 


93 (fixed) 


7.3a 




D* + n 


- (E) 


1.18 ±0.05 


71.4±1.7±7.3 


2608.7(fixed) 


93 (fixed) 




D(2750)° 


D*+n 


" (E) 


1.23 ±0.07 


23.5±2.1±5.2 


2752.4±1.7±2.7 


71±6±11 


4.2a 


D*(2760) u 


D+ty- 


(A) 


1.41 ±0.09 


11.3±0.8±1.0 


2763.3±2.3±2.3 


60.9±5.1±3.6 


8.9a 


D3(2460) + 


D%+ 


(B) 




110.8±1.3±7.5 


2465.4±0.2±1.1 


50.5(fixed) 




D*(2600) + 


D%+ 


(B) 




13.0±1.3±4.5 


2621.3±3.7±4.2 


93 (fixed) 


2.8a 


D*(2760) + 


D%+ 


(B) 




5.7±0.7±1.5 


2769.7±3.8±1.5 


60.9(fixed) 


3.5a 



fixing the parameters of all other BW components to the 
values determined in the previous fits. This final fit (Fit 
E), shown in Fig. [3] (bottom), has a x 2 / NDF of 244/207. 

Systematic uncertainties on all fit results are estimated 
by varying the parameters that were fixed in the fits 
and by varying the bin width and mass range of the his- 
tograms. In addition, the BW shape used for the new sig- 
nals is replaced by that for a D-wave decay, and we vary 
the background model according to deviations observed 
when this model is used to fit the smooth distribution in 
the wrong-sign samples. A systematic uncertainty is also 
estimated from a possible contribution of the D' 1 (2430). 
Finally, we estimate uncertainties on the mass values due 
to uncertainties in the magnetic field and the SVT mate- 
rial density. Effects due to possible interference between 
the decay amplitudes for different excited states and the 
background amplitudes are ignored in this inclusive anal- 
ysis. 

The final model for the M(D* + ir~) distribution is used 
to extract the signal yields as a function of cos 9h ■ We di- 
vide the data into 10 sub-samples corresponding to cos Oh 



TABLE II: Properties of the predicted states [J. The value 
of the parameter h depends on the state. 



State 


Predicted Mass 


J y 


cos8h Distribution 


D&(2S) 


2.58 GeV/c 2 
2.64 GeV/c 2 
2.44 GeV/c 2 
2.40 GeV/c 2 

2.49 GeV/c 2 

2.50 GeV/c 2 


0" 


oc cos 2 9h 


Df(2S) 
Di(lP) 


1" 


oc sin 2 Oh 


1+ 


oc 1 + h cos 2 Oh 


dUip) 


0+ 


decay not allowed 


Df(lP) 


1+ 


tx 1 ± h cos 2 9h 


DlO-P) 
Dl{lD) 


2+ 


oc sin 2 6h 


~2.83 GeV/c 2 
2.82 GeV/c 2 


2" 


oc 1 ± h cos 2 Oh 


D\{1D) 


1" 


oc sin 2 Oh 


dUid) 


~2.83 GeV/c 2 
2.83 GeV/c 2 


2" 


oc 1 ± h cos 2 Oh 


d 3 3 (id) 


3" 


oc sin 2 Oh 



intervals of 0.2 between —1 and ±1. Each sample is fitted 
with all shape parameters fixed to the values determined 
above. The yields extracted from these fits arc plotted 
for each signal in Fig. 01 For the i?i(2420) we measure 
the helicity parameter h = 5.72±0.25, where the error in- 
cludes both statistical and systematic uncertainties. This 
value is consistent with the measurement by ZEUS [Hj]. 
The cos8 H distributions of the D^(2460) and D*(2600) 
are consistent with the expectations for natural parity, 
defined by P = (— 1) J , and leading to a sin 2 9h distribu- 
tion. This observation supports the assumption that the 
enhancement assigned to the D*(2600) in the D + ir~~ and 
D* + ir~ belong to the same state; only states with nat- 
ural parity can decay to both D + ir~ and D* + n~ . The 



distribution for the D(2550)° is consistent with 
2 °h as expected for a J p = 0~ state. 

D+TT-) 



COS t?H 

pure cos 

The ratio of branching fractions ^(D' * ^d* +lr -) ( wnere 
D** labels any resonance) can be useful in the identi- 
fication of the new signals with predicted states. We 
compute this ratio for the L>|(2460) , £>*(2600)°, and 
L*(2750)° using the yields obtained from the fits to the 
total samples and correcting for the reconstruction effi- 
ciency: (Ndtt/£dtt)/{Nd'tt/£D'tt)- The efficiencies and 
yields are shown in Table Q] We find the following ratios: 



5(Z^(2460)° D+tt-) 
B(D*(2A60)° D*+tF) 
S(D*(2600)° -> D+tt-) 
B(D*(2600)° D*+n-) 
B(D*(2760)° -> D+n- 
B(D(2750)° -4- D*+tt~ 



= 1.47 ±0.03 ±0.16, 



= 0.32 ±0.02 ±0.09, 



= 0.42 ±0.05 ±0.11. 



The first uncertainty is due to the statistical uncertainty 
on the yields. The second uncertainty includes the sys- 
tematic uncertainty on the yields, the systematic uncer- 
tainty due to differences in PID and tracking efficiency, 
and the errors from the branching fractions for the de- 
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,Y°c [l+(5.72±0.25) 
xcos 2 (9 H )]£ H 



D,(2420) 



cos9„ 
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Y « cos 2 (8 H )e H 

D(2550) 



X 



cos6 u 
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in 2 (8 H )E H \ 



K 



COS0 H 




D(2750) 

[l-(0.33+0.28) 
xcos 2 (9 H )]£ H 



cos6„ 



FIG. 4: (color online) Distribution in cos 6h for each signal in D* + tt . The error bars include statistical and correlated 
systematic uncertainties. The curve is a fit using the function Y shown in the plot; eh is the efficiency as a function of cos#h. 



cay chains Q. Although in the last ratio the signal in 
the numerator may not be the same as the signal in the 
denominator, we determine the ratio, as it may help elu- 
cidate the nature of this structure. 

In summary, we have analyzed the inclusive produc- 
tion of the D + ir~, D°ir + , and D* + ir~ systems in search 
of new D meson resonances using 454 fb _1 of data col- 
lected by the BABAR experiment. We observe for the first 
time four signals, which we denote Z?(2550)°, £>*(2600)°, 
D(2750)°, and D*(2760)°. We also observe the isospin 
partners J D*(2600)+ and £>*(2760)+. The L>(2550)° and 
£>*(2600)° have mass values and cos Oh distributions 
that are consistent with the predicted radial excitations 
D\(2S) and D\(2S). The L>*(2760)° signal observed in 
D + tt~ is very close in mass to the Z?(2750)° signal ob- 
served in D* + it~] however, their mass and width values 
differ by 2.6a and 1.5a, respectively. Four L = 2 states 



are predicted to lie in this region [1], but only two are ex- 
pected to decay to D + n~ . This may explain the observed 
features. 
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